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James R. Stone* 

ABSTRACT 

A one-dimensional-flow model is developed which relates the 
pressure drop for subcooled boiling inside straight, circular 
tubes with constant heat flux to parameters generally known (phys- 
ical properties, heat flux, mass velocity, and geometry), and one 
unknown variable, the net fraction of heat added to the fluid 
which goes into vaporization. The basic equations of change are 
applied to a differential control volume. It is assumed that a 
modified single-phase friction factor may be used. Thus, an equa- 
tion for the local pressure gradient is obtained, dependent only 
on known quantities and the vaporization-rate parameter. This 
equation is then integrated for constant heat flux, physical prop- 
erties, and vaporization-rate parameter, The overall pressure 
drop is thus obtained as a function of known quantities and the 
effective mean value of the vaporization-rate parameter. The meth- 
od of correlation is to determine mean values of this vaporization- 
rate parameter from experimental pressure drop data and then cor- 
relate this parameter as a function of test variables 

The effects of non-equilibrium void fraction on pressure drop 
are adequately approximated to give a satisfactory correlation of 
the pressure-drop data. However, since the mean vaporization rate 
parameter is really a lumped pressure-drop parameter, based on sev- 
eral simplifying assumptions, this correlation should not be used 
to predict void fraction. The correlation is based on data for 
water at pressures from 17 to 400 psia, mass velocities from _ 
0.61x10° to 10.4x10® lbm/(hr) (ft^) , and heat fluxes from 0.13xl0 b 
to 3.45x10® Btu/(hr) (ft^) . Calculations based on this correlation 
are simplified since no independent prediction of void fractions or 
experimental void-fraction data are required in order to predict 
the pressure drop. 


INTRODUCTION 

This analysis deals with the pressure drop during subcooled 
boiling at low pressures in straight circular tubes with constant 
heat flux. Subcooled boiling generates vapor within a fluid whose 
bulk temperature is below its saturation temperature. Thus, when a 
cold liquid contacts a sufficiently hot surface, vapor forms at the 
hot wall, and some of it condenses in the cold liquid stream. The 
resulting increased volume and velocity of the stream produces both 
a higher heat transfer coefficient and a pressure drop often several 
times greater than that for all-liquid flow. 

Knowledge of pressure drop in subcooled boiling is important 
in the design of compact power systems with high-heat-flux boilers 
and direct-boiling nuclear reactors. The pressure drop must be 
known to determine local saturation temperatures and pumping power 
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requirements. A better understanding of the subcooled-boiling re- 
gime also contributes to the understanding of boiling dynamics. 

The problem of subcooled-boiling pressure drop is a particu- 
lar area of the more general problem of pressure drop in two-phase 
flow, which has been the subject of numerous experiments and cor- 
relations. A comprehensive review of these studies is presented 
by Tong [ 1 ] . In order to apply such two-phase pressure-drop cor- 
relations to subcooled boiling, it is necessary that the vapor 
fraction be known. Since this is not generally known, especially 
at low pressures, most existing correlations of subcooled boiling 
pressure-drop data have been made on a purely empirical basis, 
and the extent of their validity is unknown. Three of the typical 
correlations are: Reynolds [2] for pressures from 45 to 100 psia, 

Owens and Schrock [3] for 50 to 400 psia, and Mendler, et al. [4] 
for 800 to 2150 psia. Thom, et al. [5] have recently developed a 
relation for the mean density of water in subcooled boiling at 
750 to 1000 psia$ this relation and Thom's earlier two-phase 
pressure-drop correlation [6] were used successfully to calculate 
subcooled-boiling pressure drop. Although this approach is more 
basic than the preceding correlations [2, 3, 4], the results may 
be limited to high pressures. 

The differences in form of the various subcooled-boiling 
pressure-drop correlations and the fact that some data, such as 
those of Jeglic, Stone, and Gray [7], do not agree well with any 
of these correlations indicate that there is some uncertainty 
about the proper method of analyzing the data. Thus, it is desir- 
able to develop equations which will allow the correlation of data 
over a wide range of test conditions. The purpose of this study 
is to produce a broadly applicable, but minimally complicated, 
correlation of subcooled-boiling pressure drop for low-pressure 
water flowing in tubes with constant heat flux. One of the avail- 
able void-fraction predictions, e.g. [8, 9, 10], could be used, 
but this would greatly increase the complexity of the analytical 
formulation. However, the pressure-drop formulation must take 
into account the effects of void fraction, in order to give a 
satisfactory correlation of the pressure-drop data. This is done 
by determining a mean vaporization-rate parameter from the 
pressure-drop data to account for the phase non-equilibrium. Data 
are correlated for water at pressures from 17 to 400 psia, mass 
velocities from 0.61x10® to 10.4X10® lbm/(hr) (ft^) , heat fluxes 
from 0.13x10® to 3.45x10® Btu/(hr) (ft^) , tube diameters from 
0.118 to 0.375 inch, and tube length-to-diameter ratios from 25 
to 127. . 


DERIVATION OP CORRELATING EQUATIONS 

The basic equations of change for homogeneous two-phase flow 
(mean gas velocity equal to mean liquid velocity) are applied to 
the differential control volume shown in Fig. 1. The void frac- 
tion is formulated in terms of the mean local velocity. This 
velocity is then related to the heat balance by means of a 
vaporization-rate parameter r, which is the net fraction of the 
heat added through the wall that goes into vaporization. The 
resultant expressions for velocity and void fraction are substi- 



tuted in the pressure-gradient equation, giving a relation which is 
a function of known variables, y, and the two-phase friction factor., 
fipp . The pressure-gradient equation is then integrated, assuming 
heat flux, physical properties, fTP, and r are constant. Assuming 
a simple relation for the effective mean two-phase friction factor, 
fipp, the overall pressure drop is obtained as a function of known 
variables and the effective mean value of the vaporization-rate 
parameter, r. The method of correlation is to determine effective 
mean values of the vaporization-rate parameter 7 from experimen- 
tal pressure-drop data and then correlate 7 as a function of 
known parameters. 

Continuity Equation 

The continuity equation is used to obtain a relationship be- 
tween the local average void fraction and the local average veloc- 
ity. The continuity equation for steady two-phase flow may be 
written as follows s 


d(v gPg A g + v z p z A z ) 


= 0 


( 1 ) 


Introducing the assumptions of equal average phase velocities, i.e. 
Vg = V£ = v, and constant densities, and defining the average void 
fraction a = A g /A, there results: 


da Pi ~ ^(p l ~ Pg) 
d v " v( Pl - p g ] 


( 2 ) 


Noting that the density of the liquid is much greater than that of 
the gas, Eq. (2) simplifies to: 


da __ 1 - a 
dv “ v 


Separating variables and integrating, noting that at a =-<0, , 
v = G/p j , yields 


a = 1 - G/p jV 


(4) 


Vapor Generation Equation 


In this section the local mean velocity is related to the heat 
balance by means of the vaporization-rate parameter. As heat is 
added to a fluid flowing in the subcooled-boiling regime, an amount 
of fluid per unit time dW v is vaporizing at the hot wall, and an 
amount dW c is condensing in the subcooled stream; thus the net 
rate of vaporization in the control volume, dW v - dW c , is equal to 
d( P gvA ). Prom the continuity equation (Eq. (4)), the net vapori- 
zation^rate is given by 


A d(pgVa) 



(5) 


The maximum rate of vapor formation is given by the rate of heat 
input, 7 t Dq dz, divided by the enthalpy of vaporization A. The 
ratio of the actual net rate of vaporization to the maximum is 
given by the expression: 



r = 


( 6 ) 


DAp 

"4q 


g dv 
dz 


By rearranging Eq. (6), the rate of change of velocity with dis- 
tance may be obtained as follows: 


dv 4qr / 7 \ 

dz DAp g v ' 

Bowring [8] reports that the fraction of heat going into vaporiza- 
tion is essentially constant. But Levy [9] indicates that this is. 
not the case and instead proposes an exponential variation of r 
with heat-balance quality. This relationship is then shown by 
Kroeger and Zuber [10] to be an approximation. For simplicity in 
the present formulation, a mean, or lumped, value will be used. 
Thus, assuming that r is constant, Eq. (7) may be integrated, 
noting that v = G/pi at z = L 0 (the inception of boiling). 

Thus v is obtained as a function of z. 


Momentum Equation 


v 


4qr / 

l DAp ^ 


L 0 ) 


( 8 ) 


A momentum balance about the differential volume yields the 
following: 



where F is an orientation factor equal to the cosine of . the angle 
between the +z axis and the gravity vector. Again assuming con- 
stant densities, and making use of the continuity equation (Eq. (4)), 
the momentum equation may be simplified to yield the following 
equation for the pressure gradient: 


dP 

dz 



_G_ dv 
g c dz 



a) 


( 10 ) 


Two-Phase Friction Factor 


The two-phase friction factor is defined, by analogy to the 
single-phase (Fanning) friction factor, as the ratio of twice the 
shear stress to the total momentum flux or the resistive force di- 
vided by surface area times dynamic pressure: 


f 


tp 


( p 


2 t Sc 

7V 2 + p s v 2 y A 


(nr 


Noting that p g /p z Is much less than unity, and using Eq. (8), the 
frictional pressure drop term may be written 


il s 
D 


2f 

Dg 


^•(Gv + PgV 2 ) 


( 12 ) 



The pressure gradient may now be obtained in terms of known param- 
eters and r, substituting Eqs. (4), (8)* and ( 12 ) in Eq. (10), 
assuming 1 + (p /p^) = 1 and 1 + 2{pg/p l ) = 1. 



Integrated Pressure Drop Equation 


The pressure gradient equation may be integrated assuming r 
and f^p are constant (as well as q* G* D and properties). In- 
tegrating from z = L Q to z = L* and defining = L - L 0 * the 
following equation for the pressure drop is obtained. 



The bars are written over r and fp p to indicate that these are 
effective mean values. Equation (14/ may be nondlmensionalized as 
follows. The dimensionless group* q /(j~K, is commonly called the 
boiling number and is denoted in the remainder of this report by 



Evaluation of Vaporization Rate Parameter 


In order to make use of the model presented herein to calcu- 
late pressure drops* the mean vaporization-rate parameter r must 
be known for the conditions of interest. To develop a means of 



predicting r, the data of Reynolds [2] and Jeglic, Stone, and 
Gray [7] are examined. For the purpose of evaluating y from 
pressure drop data, Eq. (15) must be solved for r. To get this 
solution in closed form, a series approximation for the logarith- 
mic term is used. Although this introduces some inaccuracy, the 
gravitational pressure drop term is small compared with the iner- 
tial and frictional term, so that little error in r is intro- 
duced by this approximation. Thus, Eq. (15) may be rewritten as a 
quadratic in r and solved as follows: 


r = 


JlY 


4 . £ _ _ 2 _ 

,2a/ a 2a 


where 


a = 


'M 3 


2/ L b 


¥ f tp N b(i=rJ +! FN b^ 


P z L b (g/g c ) 
.° 2 /p g g o . 


(16) 


(16a) 


J b 


V 


4N blTTi + 4 f*tp N b\TT 




P z L b (g/g c ) 


a2 /PgSo 


(16b) 


AP t 


c = 


(GVp g g c ) 


- 2f . 


-F 

P z L b ( g /g c ) 

WW 

Cr^/ PgS c 


( 16c) 


In order to solve for r from Eq. (16), f b must be known. 
Thom [6] found that the two-phase friction factor is not greatly 
different from the liquid-phase friction factor at low vapor qual- 
ities. Therefore, a simple relation of the type, f b p « Re^O-2, is 
chosen. Several trials indicate that the following relation ap- 
plies . 

f tp = 0.020 Re~° o2 (17) 

-0 2 

Higher values of the_coef f icient multiplying Re^ yield in- 

creased scatter in r values calculated from experimental data. 

The difference of this equation from a standard single-phase cor- 
relation may well be due to the combined effects of the definitions 
and simplifying assumptions used. 


RESULTS AND DISCUSSION 


The subcooled-boiling pressure-drop data of Reynolds [2], and 
Jeglic, Stone, and Gray [7] are now examined,, The experimental 
pressure-drop data yield mean values of the vaporization-rate pa- 
rameter r from Eq. (16). As noted by Kroeger and Zuber [10], 
the point of boiling initiation is quite important, and no widely 
applicable means is available to determine L Q . Therefore, the 
only data used herein are those for which are given either L 0 [2] 
or the wall-temperature distribution, from which L 0 may be esti- 
mated [7], Whether this L indicates the first surface bubble 
nucleation or the point of the first bubble departure is uncertain. 



The mean vaporization-rate parameter is correlated as a function of 
known variables. Finally* as a check* the correlation obtained is 
used to predict the pressure drop* and this calculated pressure 
drop is then compared with experimental data. This is done for the 
data [2*7] used in obtaining this correlation and also for the data 
of Owens and Schrock [3]. 

Correlation of Mean Vaporization-Rate Parameter 


Figure 2 shows the mean vaporization-rate parameter r plot- 
ted against subcooling number N sc for various values of boiling 
number Nv over a wide range of test variable^ at an exit pres- 
sure of 100 psia [2*7]. It can be seen that r decreases with in- 
creasing N sc for constant N^ Also* r increases with increas- 
ing N^ at constant N„ c . No effects of heat flux or mass veloc- 
ity* except as accounted for in N^* are seen* nor are there any 
effects of geometry not accounted for in the equations. The data 
for a given pressure may be reduced to a single curve by plotting 
X against N SC N£°° * as shown for 100 josia data [2*7] in Fig. 3. 

To account for the effect of pressure* r is plotted against 
Nsc^b ' / (Pz,/Pg ) 0 ’ 5 in Pig- 4; data for very low pressure drops 
are not shown on this figure since there is considerable scatter in 
such data. It should be noted that N sc is based on heat-balance 
enthalpies. As has been previously observed [8*9*10]* liquid tem- 
peratures in the subcooled-boiling regime are less than a heat bal- 
ance would indicate; thus* establishment of limiting behavior can 
only be approximate. In the limit as N sc -*■ 0* r is assumed to 
approach 1 . 0 ; so the following simple equation yielding this limit 
is used: 


[_1 + 0.0145 N sc N£°- 7 (p z /p g ) 0 * 5 


( 18 ) 


No attempt is made to establish limiting behavior at low r. This 
would require predicting the inception of boiling and is beyond 
the scope of this report. 


Comparison of Experimental and Calculated Pressure Drops 

Pressure drops* calculated from Eq. (15) with f^p from 
Eq. (17) and r from Eq. (18)* are compared with experimental data 
in Fig. 5* where (APb) Ex p is plotted against (AP] 3 ) 0 a ]_ . The data^ 
used to obtain the correlation and some lower pressure-drop data 
[2,7] are shown in Figs. 5(a) and (b). Figure 5(a) shows data for 
water boiling in a 0.375-inch-diameter tube at pressures from 45 to 
100 psia with mass velocities from 1 . 56 x 106 to 2 . 34xl0 6 lbm/(1%r)(ffc 2 ) 
and heat fluxes from 0.13x10 s to 0.30x10 s Btu/(hr) (ft 2 ) [2]. Per- 
haps more importantly* the boiling number Nfo ranges from 
0.060x10“^ to 0.22x10"3. All the data agree within ±30 percent 
±0.1 psi. (The absolute value is included in the scatter band de- 
scription since the low AP-^ data are difficult to measure accu- 
rately. ) Figure 5(b) shows data for water boiling in 0.23-inch- 
diameter tubes at pressure from 17 to 100 psia with mass velocities 
from 1.04x10 s to 10.4x10 s lbm/(hr) ( ft 2 ) and heat fluxes from 
0.52x10° to 3.45x10° Btu/(hr) (f t 2 ) [7]; N^ ranges from 0.22x10"^ 
to 1.56X10 -3 . 96 Percent of the data agree within ±30 percent 

±0.5 psi. As a further check of the correlation* the data of 



Owens and Schrock [3], not used in developing the correlation, are- 
shown in Fig. 5(c) for Nk ranging from 0.17xl0“3 to 0.87xl0 _ 3. 
These data are for water 'boiling in 0. 1181-inch- and 0.1824-inch- 
diameter tubes at pressure from 50 to 400 psia with mass velocities 
from 0.61x10® to 3.9x10° lbm/fhr) (f t 2 ) and heat fluxes from 
0.21x10® to 1.27x10® Btu/(hr) (ft 2 ) . 98 Percent of the data agree 

within ±30 percent ±0.25 psi. It should be noted, for comparisons 
regarding scatter, that Reynolds [2] correlated his data within 
±20 percent and that Owens and Schrock [3] correlated their data 
within +27 to -31 percent. Therefore, considering the wider range 
of test variables, especially boiling number and pressure, the 
data scatter obtained with this model is considered acceptable. 
Thus, the non equilibrium homogeneous -flow model presented herein 
is, at least effectively, valid over the ranges of variables given 
above. 


SUMMARY OF RESULTS 

A one-dimensional, nonequilibrium, homogeneous-flow model is 
developed herein which relates the pressure drop in subcooled boil- 
ing inside tubes with constant heat flux to variables generally 
known (physical properties, mass velocity, heat flux, and geometry) 
and one unknown variable, the net fraction of the heat added that 
goes into vaporizing liquid. This variable is referred to herein 
as the vaporization-rate parameter. A correlation of the mean val- 
ue of the vaporization-rate parameter is presented. The data cor- 
related are for water at pressures from 17 to 400 psia, mass veloc- 
ities from 0.61x10® to 10.4x10® lbm/(hr) (ft 2 ) , and heat fluxes 
from 0.13x10® to 3.45x10® Btu/(hr) ( ft 2 ) y the boiling number, de- 
fined as heat flux divided by the quantity, mass velocity times la- 
tent heat of vaporization, ranges from 0.06xl0“3 to 1. 56x10“*. The 
effects of void fraction are lumped in the vaporization-rate param- 
eter. Therefore, this model should not be used to predict void 
fraction explicitly. Calculations based on this correlation are 
simplified since no independent prediction of void fractions or 
experimental void-fraction data are required in order to predict 
pressure drop. 


NOMENCLATURE 

A cross-sectional area for flow, in. 2 

a parameter defined in Eq. ( 16a) , dimensionless 

b parameter defined in Eq. (16b) , dimensionless 

c parameter defined in Eq. (16c), dimensionless 

D tube inside diameter, in. 

F orientation factor, dimensionless 

f friction factor, dimensionless 

G mass velocity, lbm/( sec) ( in. 2 ) 

g acceleration due to gravity, in. /sec 2 
g c conversion factor, 386 (ibm) ( in. )/( lbf ) ( sec 2 ) 
h enthalpy, Btu/lbm 

L length, in. 

N-j-j boiling number, q/GA, dimensionless 

N sc subcooling number, [h s i - (h 0 + h e )/2]/A, dimensionless 
P pressure, psia 

q heat flux, Btu/( sec) ( in. 2 ) 

Re-^ liquid Reynolds number, DG/pj, dimensionless 



v velocity, in. /sec 

W rate of phase change, Ihm/sec 

z axial distance along boiler tube, in. 

a void fraction, Ag/A, dimensionless 

T vaporization rate parameter, dimensionless 

A difference, value at inception of boiling minus value at tube 
exit 

A enthalpy of vaporization. Btu/lbm 

(i viscosity, lbm/( in. ) ( sec) 

p density, lb/in.3 

t shear stress at wall, psi 

Subscripts : 

b boiling 

Cal calculated 

c condensing 

Exp experimental 

e tube exit 

g gas phase 

l liquid phase 

o inception of boiling 

s value at saturation temperature 

tp two-phase 

v vaporizing 
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Figure I. - Differential control volume for subcooled boiling 
flow. 
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Figure 2. - Mean vaporization-rate 
parameter as a function of subcool- 
ing number for various boiling 
numbers; exit pressure, -100 psia. 
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Figure 3. - Mean vaporization-rate param- 
eter as a function of dimensionless pa- 
rameter N SC N^°* 7 for exit pressures of 
17 and 100 psia. 
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(A) DATA OF REYNOLDS (REF. 2). 

Figure 5. - Comparison of experimental and 
calculated pressure drops. 
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(B) DATA OF JEGLIC, STONE, AND GRAY (REF. 7). 
Figure 5. - Continued. 
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(C) DATA OF OWENS AND SCHROCK (REF. 3). 

Figure 5. - Concluded. 



